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ABSTRACT 

T— I ■ 

^ ' Context. Satellite accretion events have been invoked for mimicking the internal secular evolutionary processes of bulge growth. However, 
N-body simulations of satellite accretions have paid little attention to the evolution of bulge photometric parameters, to the processes driving 
this evolution, and to the consistency of this evolution with observations. 

Aims. We want to investigate whether satellite accretions indeed drive the growth of bulges, and whether they are consistent with global scaling 
relations of bulges and discs. 

Methods. We perform A'-body models of the accretion of satellites onto disc galaxies. A TuUy-Fisher (M oc V^J^} scaling between primary and 
satellite ensures that density ratios, critical to the outcome of the accretion, are realistic. We carry out a full structural, kinematic and dynamical 
■~~>^ ' analysis of the evolution of the bulge mass, bulge central concentration, and bulge-to-disc scaling relations. 

Results. The remnants of the accretion have bulge-disc structure. Both the bulge-to-disc ratio (B/D) and the Sersic index (n) of the remnant 
bulge increase as a result of the accretion, with moderate final bulge Sersic indices: n = 1.0 to 1.9. Bulge growth occurs no matter the fate of 
^ ' the secondary, which fully disrupts for ajp = 3 and partially survives to the remnant center for q-tf= 3.5 or 4. Global structural parameters 
i—{ ' evolve following trends similar to observations. We show that the dominant mechanism for bulge growth is the inward flow of material from 
C/J the disc to the bulge region during the satellite decay. 

Conclusions. The models confirm that the growth of the bulge out of disc material, a central ingredient of secular evolution models, may be 
^ triggered externally through satellite accretion. 
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. N-body simulations 
I 

1. Introduction of several Gyr, see KK04; lAthanassoulal l2005l) have been 

termed secular evolution processes. They include bulge growth 

Both observationally and theoretically, the central bulges of ^^^^ -^^^^^^^ ^i^^ instability followed either by a verti- 

disc galaxies have received much atte n tion since the last decade insta bility in the bar itsel f ( Hohl 1971; Combes . & Sanders, 

hm/f'^' ^ i Kormendy& KennicuttI ^ Ipfenni.er ^ Norman I behattista et al.l Boi. or 
12004. hereafter KK04) . However, we are far from understand- 



by bar dissolution due 


to a central mass concentration 


(Hasan & Norman 1990; 


Hasan et al. 


1993). KK04 hst satel- 



ing the processes that led to their formation. Several see- 

narios have been proposed. The similarities between bulges ijte accretion as a secular process as well, ttiggered by the en- 
and ellipticals have prompted models of fast, e arly forma- ^-^^^^^^^^ galaxy rather than internally. Bulge forma- 
tion fromih^.col^ halo fegenetalj tion/growth co uld come fro m the accreted sate llites themselves, 
lm2; Gilmore & Wvse 1998'), or as the result of violent in- ^^^^ ^j^^ jpfennigerlll 993t iMihos et al jfTigi, or from a 
teractions and major mergers of collap sed pre-galactic clumps .oj^bination of both dAmierri et al.l2001. hereafter ABPOll) . 
prior t o the formation of the thin disc llKauifmann et alJll993[ 

I1994I |Kaufrmann"l996'). But the realization that the bulges Minor mergers have been extensively studied using 



1993 



2005 



are old paradigm is not entirely true (Kormendv 1993) has numerical simulation^n the past j^uinn et alj 
led to the development of models in which bulges grow over Velazquez & White 1999; Bournaud et al. 2004, 
long timescales, after the formation of the disc. Processes that l Agpvps ^ Vglazque7„200.5a.b , ; Aceves et al. .2QQ ^ to name just 
take many galactic rotation periods to complete (of the order ^ f^^). However, interest in satelhte accretion has focussed 
more on the effects o n the disc, such as thickening and warpin 



more on ttie etiects o n ttie disc, sucti as ttiickening and warping 
dWal ker et all Il996t iHuang & Carlberg . .1997; iHelmi et all 
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balcells@iac.es; jalfonso@iac.es; cglez@iac.es 1999.) . with httle attention given to the evolution of the 
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bulge photometric parameters and to the processes that drive 
this evolution. Simulations can easily tell us whether minor 
mergers can indeed drive the growth of early-type bulges. 
They can also tell us whether infall-driven bulge growth is 
consistent with global scaling relations of bulges and discs. 

ABPOl investigated the previous questions with N-body 
simulations of the accretion of dense, spheroidal satellites. 
After the accretion of the satellite, bulge-disc decompositions 
showed that the bulge grew and the Sersic index n increased 
proportionally to the satellite mass. However, most ABPOl 
satellites were dense spheroids, not representative of the ma- 
jority of satellites of disc galaxies in the local Universe; those 
satellites reached the center of the remnant with little tidal strip- 
ping, yielding bulge growth measures that may have been unre- 
alistic. Accretion-driven bulge growth is sensitive to the density 
contrast between the primary and the merging satellite. While 
high-density satellites such as those of the ABPOl study de- 
posit most of their mass, orbital energy and angular momen- 
tum in the center of the remnant, low-density satellites disrupt 
during the merger, and their mass gets deposited at intermedi- 
ate radii. The fate of satellites with realistic densities cannot be 
ascertained a priori or with analytical means, but it can with 
simulations. 

We have expanded on ABPOl work by running N-body 
simulations of the accretion of satellites that have lower den- 
sities than ABPOl experiments. Our satellites have more com- 
plex inner structure than ABPOl spheroidal satellites: they 
comprise bulge, disc and halo. Density ratios between primary 
and satellite derive from using the Tully-Fisher (TF) relation in 
the scaling of the primary galaxy and the satellites. We use our 
models to study whether satellites with realistic densities drive 
bulge growth as compared with the ABPOl experiments, and 
to study what physical processes are involved in the evolution 
of the central component of the galaxies. We will show that 
the fate of the satellites depends critically on the exponent q-tf 
of the TF relation used for density scaling: varying axp over 
a realistically motivated range the fate varies from survival to 
full disruption. Our models suggest that TF-scaled accretion 
events produce an increase of both the BID and the Sersic in- 
dex n of the bulge, as in the high-density models of ABPOl, but 
the mechanisms that drive these changes are different in both 
cases. We find that TF-scaled satellite accretion processes alter 
galaxies systematicaly, giving place to remnants whose prop- 
erties are scaled between them, depending on the mass ratios 
and the orbit of the encounter. We will show that the growth 
of the bulge out of disc material, a central ingredient of secular 
evolution models, may be triggered externally through satellite 
accretion. 

We will indistinctly refe r to our models as "satellite accre- 
tions" and "minor mergers" jBertschik & Burke rt 2003 2004). 

This paper is organized as follows. Models are described 
in ^ We explain the bulge-disc decomposition procedure in 
^3.11 while results on growth of bulges are presented in ^3.21 
and ^3.31 Population mixing is analysed in ^3.41 while section 
^3.5l comments the vertical structure of the final remnants. We 
analyse phase-space evolution and model kinematics in ^3.61 
and ^3.71 respectively. Scaling relations of discs and bulges in 
the remnants and the comparison to observed correlations are 



given in ^3.81 Limitations of the current models are analysed in 
^ We discuss the connections between remnants of satellite 
accretion and pseudobulges in ^ and finally, we list the main 
conclusions in ^ 

2. Models 

The primary galaxy model w as built using the GalactlCS code 
dKuiiken & Dubinski The model comprises an expo- 

nential disc ( Shu 1969), a King bulge (Kins 1966), and a dark 
halo built as an Evans model ( iKuiiken & Dubinski 1 1994[) . The 
core radius of the King model is 0.15 (units given below), and 
the concentration parameter is 6.7. The disc extends to 5 length 
units. It has an exponential surface density profile both radi- 
ally and vertically; scale length is /iq = 1-0, and scale height is 
ZD = 0.1 (Guthrie 1992; de Griis 1998). The chosen disc ve- 
locity dispersion makes the disc warm with a Toomre stability 
parameter of 2=1.7, enough for preventing bar instabilities in 
the disc when in isolation. Masses, radii and numbers of parti- 
cles for each component are given in Table [2 Throughout the 
paper, a gravitational constant of G- \ is used, the length unit is 
the disc scale length, and the total mass of the primary galaxy 
is equal to 6.44. Our primary galaxy matches the Milky Way 
when the units of length, velocity and mass are: R - 4.5 kpc, 
V - 220 km s"\ M = 5.1 x IO^^Mq. The corresponding time 
unit is 20.5 Myr. 

The satellite galaxy comprises a dark halo, a disc 
and a central bulge, built as a scaled replica of the pri- 
mary. A physically-motivated size-mass scaling was used, 
by imposing that the primary and secondary galaxies fol- 
low the TF relation. Fo llowing [Oonzalez-Garcia & Bal cells I 
(12005. hereafter GGB 05V homologous galaxies with same 
mass-to-light ratio fulfill L oc V"'" whenever sizes R and 
masses M scale as 

RocM^-^l"^^, (1) 

so we apply scaling in eqn.^to our models. Central densities 
scale as 

p oc M^"^, (2) 

hence, for M2/M1 < 1, P2/P1 increases as ajp increases, i.e., 
satellites scaled with ajF = 4.0 have higher central densities 
relative to the primary than those with lower ajp. 

Orbital parameters for the merger experiments, satellite 
masses and half-mass radii are given in Table|2 We ran experi- 
ments with satellite luminous masses equal to 1:2, 1:3, and 1:6 
of the primary bulge mass. These fractions correspond to lumi- 
nous mass ratios between the satellite and the primary galaxy 
of 1:6, 1:9, and 1:18, respectively. 

We use a coordinate system whose origin lies at the system 
center of mass, the initial merger orbit lies in the XY plane, and 
the satellite lies on the +X axis at the start of the simulation, 
with orbital angular momentum pointing to +Z. Orientations of 
the disc spins are given by standard spherical (6, <p) coordinates. 
Subindices 1 a nd 2 denote the primary and satellite galaxies, re- 
spectively. As Khochfar & Burkert 12006.) have shown that the 
halo spin planes and the orbital planes of the mergers between 
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Number of Particles (/lO^^) 








Primary Galaxy 






NP 


Disci 


Bulgel Halol Disc2 Bulge2 


Halo2 


MbuI 




Aloark ''Bui 


ho 


ZD 


(1) 


(2) 


(3) (4) (5) (6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 


(13) 


185 


40 


10 90 10 5 


30 


0.42 


0.82 


5.20 0.195 


1.0 


0.1 



Description of columns: (1) Total particle number (/lO^). (2) Number of primary disc particles. (3) Number of primary bulge particles. (4) 
Number of primary halo particles. (5) Number of satellite disc particles. (6) Number of satellite bulge particles. (7) Number of satellite halo 
particles. (8) Primary bulge mass. (9) Primary disc mass. (10) Primary halo mass. (11) Primary bulge half-mass radius. (12) Primary disc scale 
length. (13) Primary disc scale height. 



Table 2. Orbital and scaling parameters for the merger experiments. 



Model Code 


Msa.(Lum)/MG(Bulge) MsJMo 


ctf R&m/Rg 


Lz,0 


e 


^peri 


^peri 


^1 ^full merger 




(1) 


(2) (3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 



M2TF4D 


1/2 


1/6 


4.0 


0.40 


0.4815 


0.907 


0.73 


1.981 


30 


-66 


100 


M2TF35D 


1/2 


1/6 


3.5 


0.46 


0.4815 


0.907 


0.73 


1.981 


30 


-72 


100 


M2TF3D 


1/2 


1/6 


3.0 


0.54 


0.4815 


0.907 


0.73 


1.981 


30 


-80 


100 


M3TF4D 


1/3 


1/9 


4.0 


0.33 


0.3532 


0.900 


0.79 


1.793 


30 


-70 


100 


M3TF35D 


1/3 


1/9 


3.5 


0.39 


0.3532 


0.900 


0.79 


1.793 


30 


-80 


100 


M3TF3D 


1/3 


1/9 


3.0 


0.48 


0.3532 


0.900 


0.79 


1.793 


30 


-85 


100 


M6D 


1/6 


1/18 


3.5 


0.28 


0.1963 


0.892 


0.86 


1.606 


30 


-116 


122 


M2R 


1/2 


1/6 


3.5 


0.46 


0.4815 


0.907 


0.73 


1.981 


150 


-80 


100 


M3R 


1/3 


1/9 


3.5 


0.39 


0.3532 


0.900 


0.79 


1.793 


150 


-87 


100 


M6R 


1/6 


1/18 


3.5 


0.28 


0.1963 


0.892 


0.86 


1.606 


150 


-142 


154 



Description of columns: (1) Model code. (2) Initial mass ratio between luminous satellite material and primary bulge material. (3) Initial mass 
ratio between luminous satellite material and luminous primary material. (4) Tully-Fisher index for scaling. (5) Ratio of half-mass radii for the 
luminous distributions of satellite and primary galaxy. (6) Initial orbital angular momentum. (7) Orbital eccentricity. (8) Pericenter distance. 
(9) Pericenter speed. (10) Initial angle between the orbital momentum and the primary disc spin. (11) Time of full merger in the units of the 
simulation. (12) Total time computed in the units of the simulation. 



dark matter halos in cosmological simulations do not exhibit 
any correlation, we have selected arbitrary angles that avoid a 
perfect spin-orbit coupling in our experiments and that sample 
direct and retrograde cases. They also indicate that orbits are 
independent of the progenitor mass and major merger defini- 
tions. Therefore, for each satellite mass we ran a direct orbit 
with Gi - 30° and a retrograde orbit with 0i - 150°. In all 
cases, <pi - 0°, 62 = 25°, and 02 - 90°. All our experiments 
exhibit very high eccentricities (see Table|5J, according to the 
cosmological simulations that have found that almost half of 
all the mergers betwe en dark matter halos are ex pected to have 
occurred with e ~ 1 jKhochfar & Burkertll2006 l). Initial orbits 
were elliptical with apocenter equal to twice the disc outer ra- 
dius (-15 length units) and pericenter equal to the disc scale 
length (~1 length unit). Each orbit was run using q'tf=3.5. In 
addition, for several of the orbit and mass-ratio combinations, 
we ran models with ottf - 3.0 and 4.0, in order to test the 
response of the system to variations in the density contrast be- 
tween primary and secondary. 



Computations were carried out with SunBladelOO ma- 
chines. Evolution was computed using t he TREECODE of 
L. Hernquist (see 'Hernauist"l987l 1 199(1 iHernauist & Katz I 
1989; Hernauist & Barnes 1990), kindly made available by 
the author The equations of moti on are integrated using a 
time-centered leapfrog algorithm JPress etalJ ll986l) with a 
variable time step. Gravitational force was softened using a 
spline kernel, with a constant softening legth e=0.02. Applying 
quadrupole-moment corrections with a tolerance parameter 
0=0.8, the code computes forces within 1% of those given 
by a direct summation. The primary and the satellite models 
were allowed to relax separately for about 100 bulge dynami- 
cal times prior to placing them in orbit for the merger simula- 
tions. Time to full merger ranges from -65 to -140 time units. 
We evolved all models for -4 halo crossing times beyond full 
merger to allow the final remnants to reach a near equilibrium 
state, except models with mass ratio 1:6, which were evolved 
for -1-2 halo crossing times. 



In the remainder of the paper, models are named with 

a three-component code: MotTFq'[D/R], where m gives the The evolution of our systems is typical of unequal disc galaxy 

bulge-to-satellite mass ratio, a denotes the TF exponent (e.g., mergers, extensive l y desc ribed in the literature going back to 

TF35 for cxf = 3.5), and "D" or "R" describe the orbit ("D" 'Too mre & ToomreU 19721) . We show M3TF35D as an example 

for direct and "R" for retrograde). Where TFa does not appear, in Figs. [J (face-on view) and |2 (edge-on view). Bulge, disc, 

TF35 must be understood for that mass ratio and orbit. and dark halo components are plotted separately for clarity. 
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Fig. 1. Time evolution of the luminous material of model 
M3TF35D. Face-on views centered on the primary are plot- 
ted. Time is given in the upper left-hand corner of each frame. 
The total time of the simulation (100 time units) is equivalent 
to ~2 Gyr, if the primary galaxy is scaled to match the Milky 
Way. First and second pericenter passages are at f = 28 and 
t = 48, respectively. The XY plane corresponds to the initial 
symmetry plane of the primary galaxy, except in the bottom 
snapshots where the galactic plane of the final remnant is plot- 
ted. Orbit is counterclockwise, prograde with the disc rotation. 
Only 1/15 of the total number of particles is plotted in order 
to avoid saturation in the plot. Black dots: Primary particles. 
Gray dots: Satellite particles. The bottom three snapshots for 
the luminous components have been zoomed-in. 



Fig. 2. Time evolution of the luminous material for model 
M3TF35D. Edge-on views centered on the primary are shown. 
Time is given in the upper left-hand corner of each frame. The 
Z direction is perpendicular to the initial galactic plane of the 
primary galaxy, except in the bottom snapshots where it is per- 
pendicular to the galactic plane of the final remnant. Orbit is 
clockwise, prograde with the disc rotation. Only 1/15 of the 
total number of particles is plotted in order to avoid saturation 
in the plot. See caption to Fig.^ 



Long tidal tails and a bridge connecting the two galaxies ap- 
pear shortly after the initial pericenter passage at f ~ 26. In re- 
sponse to the tidal force from the companion, the inner regions 
of the primary disc develop transient spiral patterns and tran- 
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sitory non-axisymmetric distortions. At the end of the merger, 
satellite mass is distributed on a wide radial range over the pri- 
mary disc. We quantify the radial distribution in the following 
subsections. Here we can already conclude that a TF-scaling 
between primary and secondary can lead to significant disrup- 
tion of the secondary. Some material reaches the center of the 
remnant, and will contribute to the post-merger bulge light, but 
a large fraction ends up in the disc region. The contribution 
of satellites to the build-up of the thick disc had been shown 
in the low-density accretion experiment of ABPOl, and al so in 
the cosm ological disc-galaxy evolution models of lAbadi et alJ 

The evolution of retrograde mergers is similar to direct 
cases, except that t idal tails are inhibited due to the lack of spin- 
orbit coupling ( Mihos & Hernauistl ll996'). and to the weak 
pericenter impulse, as pericenter separation is greater than the 
satelhte tidal radius (GGB05). 

3.1. Surface density profiles 

The accretion-driven evolution of bulge and disc structural pa- 
rameters was inferred following steps similar to those used to 
infer bulge and disc parameters of real galaxies. We derived 
azimuthally-averaged surface density profiles for the luminous 
matter, with all the merger remnants viewed face-on. We then 
fitted a combined SersicH-exponential function to the density 
profiles. The Sersic law has been used extensively to model the 
surface brightness p r ofiles of bulges (Pers ic 1968; Graham 
2001. hereafter GOll: iMollenhoff^ & Heidtl l200lt IPrieto et al 



200lUMacArthur et alJ2002l) . Its functional form is 



I(r) = /e ■ exp {-b„ [(r/re) 



i/« 



1]} 



(3) 



where is the bulge effective radius, 1^ is the surface density 
at Te, and n is the Sersic index. The factor b„ is a function of 
the concentration parameter «, which, in the range 1 < n < 10, 
may be approximated by b„ - 1.9992 n - 0.3271, with an error 
< 0.15% in the whole range (see Fig. 7 from GO 1). The expo - 
nential law provides a good model for discs dFreeman 1119701) . 
Its functional form is 



I(r) = Iq ■ exp (-r/ro) 



(4) 



being ro the scale length of the profile and /q the extrapolated 
central surface density. The code used for fitting the radial sur- 
face density profiles was the one described by GOl with seeing- 
convolution turned off. The code uses a Levenberg-Marquardt 
nonlinear fitting algorithm to locate the;^^^ minimum by chang- 
ing all model parameters. Points inside r - 0.15 (the King 
core radius of the bulge) were excluded from the fit. We tested 
that including ~ 5 more points internal to the King core ra- 
dius or excluding ~ 5 - 10 points external to it did not affect 
significan tly the results of the fits (see the study performed by 
lAceves et al..2006.) . The final face-on, elliptically-averaged ra- 
dial surface density profiles of the luminous matter for all the 
models are shown in Figs.lSji-k. Figure|3^ is the surface density 
profile of the luminous matter for the primary galaxy before the 
accretion. 



Residuals from the fits in magnitudes, defined as Am - 
2.5 log(sden / sdens fit), appear below each surface density pro- 
file plot in Fig.|3l Typical RMS is ~0.05 mag, a quite reason- 
able result compared to typical observational errors. Secondary 
minima within the ocean of can arise from all the differ- 
ent model parameter combinations. We checked for their ex- 
istence by changing the initial parameters by a factor of ~10 
typically. We concluded that the x' profiles had well-defined 
minima in all the cases. Final fitted parameters and bulge-to- 
disc mass ratios derived from the fits, w ith their errors (boot- 
strap method: lEfron & Tibshiraniri993HPress et ar.ll994 and 
references therein) are tabulated in Table|3] 

For all our experiments, we recover a two component sys- 
tem with the outer parts well fitted by an exponential profile 
(i.e. , a disc), while the inner parts are well represented by Sersic 
fits with n larger than 1. The initial model presented n ~ 0.9. 

3.2. Growtli of bulges 

Figure |4] shows bulge growth vectors in the plane log(n) 
vs.log(B/D). Each segment originates at the location of the 
pre-merger primary bulge, and ends at the location of the post- 
merger bulge in the log (n) vs. log (B/D) plane. B/D and n de- 
rive from the Sersic+exponential fits (Table O. We also show 
growth vectors for the dense satellites accretion models of 
ABPOl, as well as real bulges from the samples of Ide Jong I 
fl996) . re-analysed by GOl, and APB95. 

In all our models, satellite accretion leads to a simulta- 
neous increase of B/D and of the Sersic index of the bulge. 
ABPOl had found a similar trend for galaxies ingesting a dense, 
spheroidal satellite. The present models show that satellites 
with disc structure and TF-scaled densities also lead to a corre- 
lated increase in n and B/D. We will show in § l3.3l that the pro- 
cesses leading to the increase of n and B/D are different for the 
two sets of models. Here, we compare the two trends. ABPOl 
and our experiments show similar features: we find a tendency 
for more massive satellites to lead to higher B/D and n; and, 
the satellite density, parametrized by axp, has little effect on 
the bulge growth for the range of TF exponents studied here. 
Quantitatively however, there are differences. The present mod- 
els have lower n indices for the same final BjD. The trends of 
ABPOl models and ours seem to lie along sequences given by 
n - j3(B/D)". As we wanted to force the line log {n)-log(B/D) 
to pass through the point corresponding to the initial model, we 
have used a gradient-expansion algorithm to compute a non- 
linear least squares fit to the present models using the following 
expression: 



log 



n 

^ ini 



- a log 



(B/D) 



{B/D), 



(5) 



where (B/D) ini and ni„i represent the values corresponding to 
the initial primary galaxy (see Table^. Our TF-scaled models 
fitted eq.|3for a ^ 1.4 + 0.8, with x^ = 0.021. So, roughly, 
in the present TF-scaled models n increases by a factor of 
~ 2.6 when B/D doubles from n = 0.9 (although our experi- 
ments do not reach such a high increment in B/D, see Table|3}- 
However, the ABPOl remnants obeyed worstly the trend given 
by eq.|5] exhibiting = 0.090 for a = 1.7 + 0.7. This means 
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Fig. 3. Radial surface density profiles of luminous matter and Sersic+exponential simultaneous fits. Upper panels of each frame: 
Surface density distributions and performed fits. Fitted values of the Sersic index n, BjD, and the RMS in magnitudes of the fit 
are shown in each panel. Panels: {a) Initial primary model, {b-h) Prograde models after the merger is complete, {i-k) Retrograde 
models after the merger is complete. Open circles: Model measurements with error bars. Dashed lines: Sersic r'^" fitted com- 
ponent. Dotted lines: Exponential fitted component. Solid lines: Sum of the two fitted components. Lower panels of each frame: 
Residuals profiles, in magnitudes. Diamonds: Residuals of the fits in the above panel. Error bars of model measurements are also 
plotted on the residuals in magnitudes. The dashed vertical hues at r = 0.15 indicate the lower boundary of the fitted region. 



that, in ABPOl models, n increased by a factor of ~ 3.2 (from 
n - \) when BjD doubled. The gentler increase in n found in 
the present experiments is welcome, as it alleviates the para- 
dox posed by the suggestion from ABPOl models that expo- 
nential bulges are fragile to merging, given that galaxies with 
exponential bulges are very common. In the light of the present 
models, we see that the fragility of exponential bulges stemmed 
from the use of very dense satellites; satellites with lower den- 
sities, as in our TF-scaled models, lead to a gentler n increase as 
BjD grows. The maximum n in our set of models is n = 1.87, 
whereas ABPOl high-density experiments got up to « = 4.30. 

Observed bulges are distributed along a sequence of in- 
creasing n and BjD (Fig.|4}. Our results indicate that galax- 



ies in the lower-left of the diagram, (low-mass, n k I bulges), 
would move toward the upper-right (massive, high-n bulges) 
when accreting a low-mass galaxy. ABPOl ran an experiment 
that showed the cumulative effects of merging (see their model 
H3D3D). They found exactly the same increase in n and in BjD 
in two consecutive mergers with exactly the same mass (1:3) 
onto their primary galaxy. Therefore, if the effects of merg- 
ing are cumulative, as this experiment seems to indicate, satel- 
lite accretion provides an evolutionary link connecting low- 
«, small bulges to more massive, high-« bulges. Moreover, 
Khochfar & Silk (2005) found that minor mergers are a mag- 
nitude more frequent than major mergers, a result that supports 
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Description of columns: (1) Model code. (2) RMS of the fit in magnitudes. (3) Disc central surface density. (4) Disc scale length. (5) Bulge effec- 
tive surface density. (6) Bulge effective radius. (7) Bulge profile Sersic index. (8) Bulge-to-disc mass ratio derived from the Sersic -(-exponential 
fit. 
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Fig. 4. Growth vectors in the log(n)-log(B/D) plane for the present TF-scaled models. Right panel shows a zoom of the left 
panel. Each segment starts at the location of the original model and ends at the n and BjD derived from the two-component fit to 
the surface density profile after the merger. Errors from Table|3]are overplotted. Growth vectors of ABPOl high-density models 
are plotted for comparison. Distributions of log (n) vs. log(B/D) for bulges of several studies are drawn too. Symbols and colours 
are explained in the legend. 



the possible relevance of minor mergers in the configuration of 
actual bulges. 

We caution that our models start out with higher B/D than 
most real galaxies (Fig. |3. This offset mainly derives from 
computational limitations: the higher range of dynamical times 
in galaxies with smaller bulges lead to more time-consuming 
computations; and, accretion experiments onto smaller bulges 
require smaller satellites, which take longer to decay. While 
we expect the main trend of increasing n and B/D to remain 
vahd in the low-B/D regime, the details will probably change. 
Small, n < I bulges should have a weaker tidal field which 



probably allows for smaller satellites to reach the center undis- 
turbed. This might lead to a steeper n-\og(B/D) trend than the 
one derived from the present models (eq. |5}- Conversely, we 
expect shallower slopes for the n-\og{B/D) trend for more mas- 
sive bulges that start off with « > 1 ; the process of n growth via 
dissipationless accretion saturates at some value of n, proba- 
bly around n ~ 4 (ABPOl). Shallower growth vectors for more 
massive bulges would keep accretion-grown bulges within the 
locus of real bulges in the n-\og{B/D) diagram. 
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Fig. 5. Face-on radial surface density distributions of the various luminous components, for all the models after the merger Solid 
lines: Total luminous material. For the initial model, only luminous material from the primary galaxy is considered. Open circles: 
Primary bulge material. Dashed-dotted line: Primary disc material. Dotted line: Satellite bulge material. Dashed line: Satellite 
disc material. 



3.3. Why does n Increase 

We now analyze the radial mass distributions per component 
for the models, in order to clarify what drives the increase of 
n and BjD. In Fig.|5]panels (b)-(k), we plot the final surface 
density profiles for particles originally belonging to each mass 
component, for all of our remnants; panel (a) corresponds to 
the primary and the satellite scaled with axF = 3.5 models 
before merging. Three distinct processes modify the mass dis- 
tribution of the galaxies. First, satellite material gets deposited 
over a wide radial range in the remnants. Second, the distri- 
bution of primary bulge material becomes more extended in its 
outer parts. And, third, the primary disc material becomes more 
centrally concentrated: in all of the models, surface density pro- 
files are not exponential but curve up inwards. We show here 
that primary disc redistribution is the dominant contributor to 
the change in the galaxy profile. 



We already noted earlier in §|3]that satellite material gets 
deposited over a wide radial range. Figure|5lshows that, in the 
cktf - 3.5,4.0 cases, satellite material does reach the center of 
the remnant, while low-density (axF = 3.0) or low-mass (M6) 
satellites entirely disrupt during the merger; their mass gets de- 
posited into an extended, torus-like structure. But the contri- 
bution of this material to the total profile is small. In no case 
does the satellite material dominate the central surface density. 
Rather, its contribution to the total surface density is between 
~ 2 - 6 mag below that of the primary. Only for a^-p - 4.0, 
the r = contribution of the satelhte bulge is similar to that 
of the primary disc. To highlight the small contribution of ac- 
creted material, the surface density profile from primary mate- 
rial alone is also shown in Figure|5](solid gray lines), together 
with the total luminous profile (solid black lines). The differ- 
ence between both profiles, which measures the contribution 
of the satellite material, never reaches 0.2 mag arcsec"^. Also, 
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excluding the secondary material modifies the parameters from 
the Sersic+exponential fits by less than ~ 20%. The situation is, 
therefore, quite unlike that found for high-density, spheroidal 
satellites by ABPOl, in which the central density was domi- 
nated by accreted material. 

To investigate which of the two primary-galaxy luminous 
components has the dominant contribution to the surface den- 
sity profile change, we analyze fractional mass distributions 
(FMD), i.e., the radii enclosing given percentages of the mass. 
Figure |6l shows the FMD for the primary bulge particles (top 
panel) and for the primary disc particles (lower panel), for 
all of the models. For the primary bulge, expansion occurs in 
the outer 20% of the mass, at radii from 0.2 to ~1.0, and is 
strongest for high-density (TF4) satellites; the inner 70% of the 
mass remains almost undisturbed in all cases. For the primary 
disc, FMD changes at all radii. The inner 20% contract (radii 
r < 0.8) while the outer mass fractions tend to expand. 

Hence, the changes in the final surface brightness profiles 
are mostly driven by the redistribution of the disc material. 
By increasing its central concentration, the disc adds material 
to the bulge region and increases its central surface density; 
and, by becoming brighter at radii where the profiles change 
from being bulge-dominated to disc -dominated, the disc ma- 
terial gives the final profile a smooth curvature that leads to 
higher-n Sersic fits. The changes in the primary bulge FMD 
have a secondary effect on the final profiles, as expansion oc- 
curs at radii where their contribution to the total profile is sev- 
eral magnitudes below that of the primary disc. The physical 
process driving the expansion of the primary bulge outer re- 
gions is the absorption by the bulge material of the orbital en- 
ergy and angular momentum. Bulge expansion is greatest for 
the denser (TF4) satellites (Fig.|6|l; but, in general, it is lower 
here than in the ABPOl dense spheroidal models, due to the 
lower mass, lower density and progressive disruption of the 
satellites we have modelled. For the disc, implosion of the in- 
ner 20% of its mass is due to the effects of the transitory non- 
axisymmetric distortions dr iven by the tidal field of th e satel- 
lite, as already described bv Mihos & Hernquistllll994l) . and to 
the deposition of dark matter from the satellite halo in the cen- 
ter of the remnant. Notice also that the primary disc outer parts 
expand, contributing to the increase of the disc scale lengths 
in the final remnants (see lines for FMD > 30% in the bottom 
panel of Fig.|6|l. 

3.4. Population mixing 

We now quantify the amount of population mixing in the re- 
gion of the bulge induced by the satellite infall. For the pur- 
pose of this analysis we simply define the region of the bulge 
as r < 0.5 (cf. Fig.|3, and study the evolution of the parameter 



rii(t) = 



Miit) 



(6) 



Mu,m(t = 0) 

where M denotes mass within r - 0.5, t is time, and / = 
B 1 ,D 1 ,B2,D2 for primary bulge, primary disc, secondary bulge 
and secondary disc material, respectively. We have called 
M-Lumit - 0) the total primary luminous mass that initially ex- 
isted at r < 0.5. Therefore, T/totaiCO) = 1, but ?7totai(?finai) 1- The 



disc contribution to the region of the bulge, initially 7;di(0) - 
0.16, rises in all of the models, in proportion to the satelUte 
mass: //Di(ffinai) - 0.18, 0.19, and 0.20 for mass ratios 1:6, 
1:3, and 1:2, respectively. These numbers are slightly lower 
for retrograde orbits, and do not depend on the TF exponent. 
The bulge contribution, initially /7bi(0) = 0.84, decreases in 
all cases in proportion to the satellite mass and also to the TF 
exponent, down to 77Bi(ffinai) ~ 0.79 for the most massive and 
densest case (M2TF4D). Finally, the satellite bulge particles 
accumulate in the bulge region, whereas the satellite disc ma- 
terial is deposited mostly in regions external to the bulge. As 
expected, deposition in the bulge region is stronger for more 
massive or denser satellites, and for direct orbits. Highest val- 
ues (M2TF4D) are 77B2(ffinai) = 0.09 and 7/D2(rfinai) = 0.03. 

3.5. Vertical structure of the remnants 

Simulations have shown that merging satellites with masses 
as low as a few percent of the disc mass heat discs 
efficiently dOuinn et all Il993t iMihos & Hernaui stl 1 1994 
'Huan g & Carlberg|ll997l) . As pointed out by IToth & OstrikeJ 
( 1992), disc thickening effects put constraints on how impor- 
tant merging processes may have been in the history of disc 
galaxies. Our models are not an exception: the final disc of all 
the remnants exhibits higher vertical scale than the initial disc. 
In this subsection we provide measurements of disc thickening 
for our models. These measurements serve two purposes. First, 
extreme thickening would falsify the hypothesis that satellite 
merging is relevant for the evolution of bulges of disc galaxies. 
And, second, our analysis will highlight that accreted satel- 
htes are im portant contributors to the building of thick discs 
jAbadi et a l. 2003b). 

In Figure0we plot the medians of the final vertical distribu- 
tions (klmedian) foi" all the modcls, computcd at four fixed radial 
positions (r = 0.3,0.5, 1.0, and 3.0), and separately for parti- 
cles initially belonging to the primary disc, the satellite bulge, 
and the satellite disc. We have also plotted the scale heights 
obtained through exponential fits to the vertical distribution of 
the total luminous material at each fixed radius. We see that 
scale heights increase for all the models and all the mass com- 
ponents. Focussing on the r - I panel, for the primary disc 
particles, kLedian increases by factors of 2 to 4 from an initial 

value \z(t = 0)Uedian ~ 0.05. 

For the satellite particles, Fig.0 shows that Izlmedian at r > 
0.5 is higher that klmedian for primary disc particles. I.e., the 
material accreted from the satellite disc contributes to build- 
ing a thick disc. In contrast, the satellite bulge particles are 
confined for all the remnants in a central (R < 2), thin (low 
klmedian) distribution. We will show in ^3.6l that such structures 
are dynamically-cold inner discs. 

Figure also shows that, for particles from the primary 
and secondary discs, scale heights increase with radius (note 
that the ranges plotted in the vertical axes change from panel 
to panel). I.e., our final di scs are flared. Disc flares a nd warps 
are common in real discs dde Griis & Peletie'r1ll997l) . and can 
be easily reproduc ed in simulations of merging small satellites 
dOuinn et alJl993h . 
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Fig. 6. Radii enclosing a given % of the mass of the primary galaxy for all the models. The abcissa values are the model codes 
from Table 12 The horizontal dotted lines are the radii initially enclosing a given % of the mass, which appears next to the line, 
on the right of the figure. Upper panel: For the particles initially belonging to the primary bulge. Lower panel: For the particles 
initially belonging to the primary disc. 



The values plotted in Fig.0 are raw measurements, which 
need to be corrected for numerical thickening which arises 
from the residual collisionality of the code. To estimate this 
correction, we assume that energy injection into vertical mo- 
tions in the disc, coming from two-body effects and from 
merger dynamics, are largely independent from each other: 
two-body disc vertical heating is dominated by interactions 
with primary halo particles, whose effects should not vary 
much if the disc is been thickened by the satellite accretion; 
and, accretion-driven injection of vertical kinetic energy onto 
the disc should not vary much under small variations of the 
disc thickness due to halo-disc two-body effects. Under that as- 
sumption, we have 

^z,.simul - -^^z.initial + ^K^,2R + ^K^,magmg (7) 

or 

<simul = ^initial + ^^l.lYi + ^^^'merging (8) 

where K, denotes vertical kinetic energy, cr- denotes verti- 
cal velocity dispersion, and initial, 2B, merging, simul de- 
note initial, two-body, merger, and final system, respectively. 
Because vertically the discs are appro ximately isothermal, h, - 
o-J{2nGpo)-^l^ tSpitzer..l94Z:.van der Kruit & Searlelll98ll) . 
where po is the in-plane volume density. In the limit of low 



mass accretion valid in the present experiments, vertical disc 
surface density is conserved, i.e., yu oc p^h- - const., hence 
/ij oc crj, and equation |S]becomes /z^.simui = /^z.initiai + ^h^,2B + 
A/zj merging- Under these approximations thickening effects add 
up linearly, i.e., 

A/Zz,merging ~ ^h~_simu\ " A/lj.,2B (9) 

where A/i,, simui is given by Fig.0 We measured A/i, 2b by evolv- 
ing the initial primary galaxy model in isolation for 100 time 
units, the typical duration of t he experiments. We obtaine d 
A/!z,2B = 0.05, a similar result to lKuiiken & Dubinski T ill 99.4 . 
From inspection of Fig. we conclude that A/zj- merging ~ 
for the low-mass satellites, where the measured thicken- 
ing is dominated by two-body effects. For most models, 

A/J2,merging//Zz,initial ~ 3, and A/z^.tnerging/Zi^anitial ~ 6 for the mOSt 

massive, retrograde experiments. 

Clearly the post-merger systems do not match thin discs 
of present-day spirals. Their properties resemble those of thick 
discs. Assuming the rebuilding of thin discs out of left-over 
gas after the merger, disc thickening is not a problem for the 
merger origin of bulge growth. The build-up of a thick disc 
may be a useful feature of our models for evolution toward ear- 
lier Hubble types, as thick discs appear to occur more often 
in early types than in later types dBurstein 111979ft . Finally, our 
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Fig. 7. Medians of the vertical distributions in the final remnants for the particles initially belonging to different components, at 
four fixed radii in the galactic plane (r = 0.3, 0.5, 1 .0, 3.0). The abcissa values are the model codes from Table|2] The horizontal 
dashed lines in each panel indicate the initial vertical scale of the primary disc at each radius. Only data points with more than 50 
particles per bin are plotted. Asterisks: Vertical scales for the particles initially belonging to the primary disc. Diamonds: Vertical 
scales for the particles initially belonging to the satellite disc. Filled circles: Vertical scales for the particles initially belonging to 
the satellite bulge. Open triangles: Vertical scales for all the luminous material. Filled triangles: Scale heights obtained through 
exponential fits to the vertical mass distribution of the total luminous matter in each remnant, at the fixed radius. 



results agree with lAbadi et alJ (l2003bl) that accreted satellites 
contribute to building the thick disc; we note, however, that the 
set of orbit inclinations covered by our experiments is rather 
limited. 

3.6. Phase-space evolution 

The structural transformations discussed in § 13.31 must be 
the result of changes in phase-space structure of the galax- 
ies. We show that this is the case by looking at Lindblad 
(energy-angularmomentum) projections of phase space. Figure 
ISlpresents specific angular momentum vs. energy diagrams (jz 
- s) before and after the accretion, for each galaxian compo- 
nent of models M3TF3D and M3TF4D. These two models 
are useful for showing the effects of satellite density (higher 
for ajp = 4.0) on the final dynamical structure of the rem- 
nants. Top panels of the figure show the initial particle dis- 
tributions of the luminous components of the primary galaxy 
and the satellite. The four intermediate panels present the fi- 
nal post-merger particle distributions of model M3TF3D, and 



the bottom panels, the same but for M3TF4D. Particle densi- 
ties in the Lindblad diagrams are depicted with logarithmically- 
spaced levels to highlight low-level structures. The z direction 
is defined by the total spin of all the luminous material of the 
primary galaxy (panels [a]-[b]); the satellite (panels [c]-[d]); 
and the final remnant (panels [e]-[h] and [i]-[l]). The thick 
dashed line in each panel indicates the location of circular 
orbits, jcirc(e), corresponding to the potential of the primary 
galaxy (panels [a]-[b]), the satellite (panels [c]-[d]), or the final 
remnant of each model (panels [e]-[h] and [i]-[l]). 

FigureOshows that material from the primary galaxy mod- 
ifies its distribution onto the s - plane only moderately. Disc 
particles, which were initially distributed in nearly circular or- 
bits, scatter over the region of direct orbits allowed by the po- 
tential (compare panel [a] with panels [e] and [i] in Fig.|8}. 
Bulge particles, with no initial rotation, acquire rotation in their 
outer regions (e > -3.0; see panels [b] and [f],[j]). Both pro- 
cesses bring the two distributions to look more like each other, 
i.e., the accretion leads to a mixing of bulge and disc in phase 
space. The changes in surface brightness profile and in the ver- 
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Fig. 8. Initial and final particle density distributions in the plane jz-s for each galactic component of models M3TF3D and 
M3TF4D, being the z-component of the specific angular momentum, and e, the specific energy. Darker grey levels in the 
density contours point to higher number of particles populating the corresponding j^-s bin (see the corresponding legend on the 
right of each panel). Contours are logarithmically-spaced. Panels (a-d): Initial material distributions in the e - jz plane of the 
two luminous components of the primary galaxy and the satellite for model M3TF3D (initial distributions for all the models are 
analogous). Panels (e-h): Post-merger distributions for the two luminous components of the primary galaxy and the satellite in 
the remnant of model M3TF3D. Panels (i-l): Post-merger distributions for the two luminous components of the primary galaxy 
and the satellite in the remnant of model M3TF4D. Dashed lines: Locus of circular orbits co-rotating (jz >0) and counter-rotating 
(jz <0), ;circ(e), in each case. 



tical structure of the galaxy (§§ 13.21 and l3.4t are consequences 
of such mixing in phase space. The transformation is gentle, as 
the global energy and angular momentum ranges do not change 
during the accretion. 

Phase-space evolution is much stronger for particles from 
the satellite. Material from the satellite disc, originally lying 
along the locus of circular orbits, end up nearly filling the avail- 
able prograde phase space, and populating some retrograde or- 
bits (compare panel [c] with panels [g] and [k]). Thus, the sec- 
ondary disc is destroyed. That orbits in the final potential are 
not circular is a result of an early disruption of the outer parts 
of the satellite. Those satellite disc particles that have not yet 
fallen onto the final remnant are laying on the tidal tails, cor- 



responding to clumps of e in the figures. Material from the 
satellite bulge shows an inverse trend: initially in a hot, non- 
rotating distribution (panel [d]) they end up clustered along 
the locus of (direct) circular orbits (panels [h] and [1]). We 
saw in § 13.51 that bulge particles lie in a flat structure after the 
merger We see now that this structure is a dynamically cold 
disc. That satellite bulge particles are deposited into circular 
orbits in the potential of the remnant indicates that dynami- 
cal friction has circularized the merger orbit prior to disrup- 
tion of the satellite bulge. Such a disc is too faint for detection 
if it exists in a real galaxy, except perhaps for inclined views. 
However, it may be brighter for other mass ratios, merger con- 
figurations, or if gas was present in the satellite that later forms 
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Fig. 9. Line-of-sight velocities of the models, using a slit placed onto the disc major axis of the remnant, for the inclined view: 
6 - 60°, 4> = 90°. The abcisas axis has been inverted in the retrograde models. The curve labelled as "satellite" comprises all the 
luminous material initially belonging to the satellite, (a) Primary and secondary initial models, {b-h) Prograde models after the 
merger is complete, (i-k) Retrograde models after the merger is complete. Asterisks: Primary disc material. Diamonds: Primary 
bulge material. Crosses: Satellite luminous material (bulge+disc). Solid line: Total luminous material (satellite+primary galaxy), 
except for panel (a), where only total luminosity material from the primary galaxy is plotted. Dashed line: Total dark matter halo 
(satellite+primary galaxy), except for panel (a), where only dark matter from the primary galaxy is plotted. 



stars jMih os & H ernqui'st1ll994l) . yielding an inner or nuclear 
disc not u nlike those found in early-type galaxies (Erwin et al. 

Our range of TF exponents yields both 
inner rings (TF=3: Fig.|S] panels [g] and [h]) and inner discs 
(TF=4: Fig.|8] panels [k] and [1]), as a result of the earlier dis- 
ruption of lower-density satellites. 

3.7. Rotation curves of the merger remnants 

Figure |9] presents the initial and post-merger rotation curves 
of the luminous components for all the experiments. These are 
means of the line-of-sight velocities measured in spatial bins 
along a virtual slit placed onto the disc major axis. We have 
used an inclined view with / - 60° {0 - 60°, (f) = 90°) instead 
an edge-on view to enhance the contribution of bulge particles 



in the center of the rotation curve. The slit length and width are 
6.0 and 0.5, respectively. 

The final rotation curves of the total luminous material 
(solid lines) are quite symmetric, with an overall shape very 
similar to the rotation curves of the early-type spirals. They 
all are very similar to the original rotation curve of the primary 
galaxy; i.e., TF-scaled satellite accretions imprint little changes 
to the original rotation curve. In the external regions of panels 
(h) to (k), rotat ion decreases outwards, as commonly observed 
in early-types ( Casertano & van Gorkom 19 91). In contrast, in 
the internal regions of the remnants, the primary bulge particles 
acquire net rotation, making the inner rise steeper than initially. 
The effect is stronger for denser satellites (compare panels [b] 
to [d] or [e] to [g]). 
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Fig. 10. Line-of-sight velocities of the particles originally in the satellite, after the merger is complete, for the inclined view: 
6 = 60°, <p - 90°. Asterisks: Satellite disc material. Diamonds: Satellite bulge material. Solid line: Total luminous material 
(satellite+primary galaxy), except for panel (a), where only total luminosity material from the primary galaxy is plotted. Dashed 
line: Primary galaxy dark material. Dotted line: Satellite dark material. 



That the inner rotation curve of the primary disc remains 
largely unmodified during the accretion is an indication that 
the merger relaxation is gentle in the present models in which 
satellite densities are set via TF. The accretion of denser satel- 
lites modelled by ABPOl violently steered the inner galaxy. 
The inner disc net rotation diminished as random motions in- 
creased, while the primary bulge, which eventually suffered a 
similar-mass merger with the satellite, absorbed much of the re- 
maining orbital angular momentum and ended up rotating even 
faster than the primary disc (cf. Fig. 8 of ABPOl). In contrast, 
the TF-scaled satellites modelled in the present paper leave the 
inner disc rotating rapidly, and imprint a weak rotation to the 
primary bulge. This rotation is faster for denser satellites of a 
given mass, which deliver more mass and angular momentum 
to the remnant center 

The details of the inner rotation curves are affected by the 
lack of rotation in the bulges of the initial primary models: 
the final inner rotation curves would be steeper if the initial 



bulges had net rotation. Similarly, the mild counterrotation of 
the bulges seen in some of the retrograde mergers (e.g., M2R, 
M3R, see Fig.|9}, would disappear if the initial bulges had net 
rotation. 

Figure [TO] presents the contribution to the final rotation 
curves of the two luminous components of the satellite sepa- 
rately. The satellite material does not contribute considerably 
to the final rotation curve in any experiment, due to its low 
surface brightness. Therefore, although the satellite bulge and 
disc material counterrotate in all of the retrograde experiments, 
the central counterrotation can not be detected in the total rota- 
tion curves of our retrograde models. This alleviates the prob- 
lem of ABPOl high-density models, which showed counter- 
rotation in excess with respect to what it i s actually observed 
dSilchenko et alJl997tlpizzella et all200l . 

Dark matter kinematics are also shown in Figs.l9l and 1101 
Figure |5] presents the rotation curves for the final dark mat- 
ter halo, while Fig. [TO] shows the rotation for the particles be- 
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longing to the primary halo and to the satelUte halo, separately. 
Analyzing dark matter kinematics is beyond the scope of this 
paper We only note that the global dark matter rotation curve 
exhibits rotation (Fig.|5Jl only in the inner parts of the mod- 
els involving more massive satellites. Initial halos were non- 
rotating. Primary dark matter acquires no rotation (Fig.llO>. 
hence the rotation feature seen in the global dark matter ro- 
tation is entirely due to (satellite) mass deposition. 

We have tested if our remnants follow the TF rela- 
tion of the initial satellites and primary galaxy. Asymmetric 
drift corrections have been applied to the maximum rota- 
tional velocities of each remnant viewed edge-on, according 
to|Binnev & Tremaine ( 1987). They are all offset to lower lu- 
minosities than those which should be expected for their ajp 
and Vcirc- Direct models exhibit differences of ~ 0.5 mag re- 
spect to the TF expected value, while retrograde cases present 
the highest dispersions (~ 1.0 mag). Nevertheless, these val- 
ues are similar to the observational dispe rsions of the TF re- 
lation in field spirals (Sakai et al. 2000. g -nvvj ~0.38 magh an d 
late-type field spirals (Matthew s et alJll998. crv = 0.82 mad). 
and field a nd clustered S O's (Neis tein et alJl999. (Ti ~0.7 mast 
[Hinz et aI1 l2001. cti ~0 .5 mag). Thus, low-satellite accretions 
could be one of the sources for the dispersion observed in the 
TF relation. 

3.8. Scaling relations of discs and bulges 

We now analyze the changes in bulge and disc global parameter 
relations in our models. Any dynamical and structural transfor- 
mation of galaxies must be compatible with the existence of ob- 
servational scaling relations. Hence, we enquire, first, whether 
accretion-driven galaxy transformations show any systematics, 
and, second, whether accretion events tend to broaden any pre- 
existing scaling relations. We will show here that many, though 
not all, of the observed scaling relations are compatible with 
the trends seen in our models. 

We use the photometric parameters in Table |3 In order to 
compare with the observed trends, the increments in the photo- 
metric parameters of the remnants are plotted vs. the increment 
in the bulge magnitude in Fig . lTTI Assuming the same M/L for 
all the bulges jPortinari et al.l2004l) . we have defined the incre- 
ment in bulge magnitude as 

AMbui = Mbui - Mbui,o = -2.5 ■ log (VVlBui/AtBui.o) (10) 

where AIbuI and AIbuI.o are the masses of the final and initial 
bulge, respectively. 

Most of the photometric parameters in our models show 
strong trends with the bulge magnitude increment (Fig.ll 1>. 
The majority increase towards higher AIbuI, e-g-, the Sersic 
index (n), the bulge-to-disc ratio (log(B/D)), the central and 
effective bulge surface densities (/io.Bui and //e.Bui), the disc 
scale lenght (/loisc), the bulge-to-disc ratio of the central in- 
tensities (log(/o3ui/^o,Di,sc)), and the central velocity dispersion 
((To); whereas the bulge eff'ective radius (log(re)), the disc cen- 
tral surface density (//o.Disc), the maximum rotational velocity 
of the disc (log (Vrot)), and the ratio between the bulge effec- 
tive radius and the disc scale length (re/ZiDisc) decrease towards 
higher AIbuI- Final discs have higher scale lengths and scale 



heights than the original primary disc, except in two of the ret- 
rograde cases (Fig.^2 panels [d] and [j]). Notice that the cor- 
responding relations for the ABPOl dense models (diamonds 
in Fig. ll 11 1 are weaker in all the cases. 

Orthogonal regressions to the distributions are shown in 
each panel of Fig. ^2 The trends exhibited by jUq.BuI, n, 
log (/(),Bui//o,Disc), log(B/D) and log(o-o) (panels [(c], [f], [g], 
[i] and [k] in Fig.ll 1> have absolute Pearson coefficients above 
0.88. The null correlation hypothesis is rejected for these fits, 
with a significance of 95% and for a two-tailed distribution; 
whereas the rest of correlations are not statistically significant. 
Nevertheless, the fits involving //g and the ratio rg/ZiDisc exhibit 
Pearson coefficients above 0.78. The correlations of Fig.^2 to- 
gether with the correlation found between log (n) and log (B/D) 
shown in eq.|5l indicate that TF-scaled satellite accretion does 
not alter galaxies randomly: changes in the bulge and disc scal- 
ing relations driven by accretion scale systematically with the 
change in bulge mass. 

Many of the changes of photometric parameters from our 
models follow trends with the bulge magnitude increment sim- 
ilar to the correlations with the bulge magnitude observed by 
BGP04b. In particular, the disc parameters (//o.Disc and /iDisc), 
central velocity dispersions, bulge central surface brightnesses, 
Sersic n and B/D ratios (compare in Fig.^Jour correlations - 
solid lines- with the slopes of the significant observational cor- 
relations from BGP04b -arrows-). 

The trends show some discrepancies with the observations 
as well. The bulge effective surface density (Fig.ll lb) and the 
ratio of the central intensities (Fig. ll lb ) correlate weakly with 
the bulge magnitude in observations (BGP04b), whereas they 
do strongly in our models (they exhibit Pearson coefficients 
> 0.8). Most relevantly, we fail to reproduce the trends in- 
volving the bulge effective radius re,Bui (Fig.ll Ih.h). rgBui de- 
creases with the accretion events, which is at odds with the fact 
that t he observed rnore luminous bulges have larger effective 
radii tHubhlell924lBinggeli et alJl "984l lMollenhoff-& Heidtl 
2001 , BGP04b). Such discrepancy diminishes if a higher mass- 
to-light ratio is assumed for the particles initially belonging to 
the bulge, while the other correlations do not change signifi- 
cantly. 

4. Limitations of the present models 

An inherent ingredient of computational work is the strong de- 
pendence of the results on the initial conditions. Our models are 
no exception, especially given the small coverage of parameter 
space. Here we discuss how the specific choices made for our 
models may have affected the results. 

1. Using the same initial primary B/D in all the experiments. 
Different initial B/D ratios would likely affect the merger 
transformation. Larger initial B/D would lead to more 
effective satellite disruption and lower changes in B/D. 
Conversely, a smaller initial B/D would lead to less effi- 
cient disruption of satellites, to stronger non-axisymmetric 
distortions and more effective inner mixing, which, to- 
gether, would have led to a faster fractional bulge growth. 
This implies that the slope determined for the evolution 
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Fig. 11. The dependence of bulge and disc parameters, and B/D ratio, on the increment of the bulge magnitude (see definition 
in the text). All parameters are derived from Table |3] Legend is the same as in Fig.0] On the left column, we have plotted 
photometric parameters obtained directly from the bulge-disc decompositions (scale lengths and central surface densities of 
bulge and disc, and the Sersic index n). On the right column, ratios between two of them, disc scale heights and typical velocities 
are shown. Solid lines: Orthogonal regressions to the model points. At the right of the identification letter of each panel, "S" 
indicates that the correlation is significant at a 95% confidence level, whereas "NS" indicates it is not. Dashed lines: Growth 
vectors of our TF-scaled models. Dotted lines: Growth vectors of the high-density models from ABPOl. Arrows: Slopes of the 
corresponding observational correlations found by BGP04b. Panels: (a) Effective radius of the (Sersic) bulge component, (b) 
Effective surface density of the Sersic component, (c) Central surface density of the Sersic component, extrapolated from the fit. 
(d) Disc scale length, (e) Extrapolated disc central surface density, (f) Sersic index n. (g) Logarithm of the bulge-to-disc ratio of 
the central intensities, (/i) Ratio between the bulge eff'ective radius and the disc scale length. (/) Bulge-to-disc luminosity ratio. (/') 
Ratio between the disc scale height at r = /id, obtained through an exponential fit, and the disc scale length, (k) Central velocity 
dispersion. (I) Maximum rotational velocity of the disc. 



in the n vs. B/D plane (eqn.|5} is valid only for the pri- 
maries with B/D similar to those of our models. Probably, 
the slope would be higher/lower for smaller/bigger initial 
B/D. But it is also probable that the slope of the growth in 
the n-B/D plane decreases as the n of the initial primary 
bulge is higher. This would make the growth mechanism 
to saturate at n ~ 3-4, preventing the departure of the 
bulges from the observational trend after several consecu- 
tive minor mergers. Understanding the behaviour for low- 



B/D is important as most early- to intermediate-type galax- 
ies have lower B/D than our models (see Fig.|4}. One par- 
ticularly interesting case is that of bulge-less primary galax- 
ies. Our models suggest that a bulge would easily develop 
in bulgeless galaxies due to the accretion of even tiny satel- 
lites, a manifestation of the well-known problem that bul- 
geless galaxies pose to hierarchical galaxy formation mod- 
els ( D'Onghia & Burkert 2004). The problems for conduct- 
ing experiments with lower B/D are essentially numerical. 
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Accretion experiments onto smaller bulges require smaller 
satellites, which imply longer orbital integrations, due to 
weaker dynamical friction, and smaller time steps, due to 
shorter dynamical times for bulge and satellite. 

2. The halo profile. 

Our haloes have been built with an isothermal profile. 
Using different initial halo profil gs,^.g. the cosrnologically- 
motivated NFW profile dNavarro et al]ll995l flQQ?!) . might 
modify the mass deposition, the stripping of the satellite, 
and the circularization of the orbits. The way in which dif- 
ferent halo profiles influence the mass deposition in the 
remnant center could provide some interesting clues about 
the central structure of dark haloes. 

3. The bulge and halo rotation. 

We used non-rotating bulges and halos in the initial mod- 
els, in order to better constrain the onset of bulge rota- 
tion due to the merger. However, r eal bulges do rotate 
(e.g. iKormendv & Illingworthiri982h . Halos probably ro- 
tate too -although httle information exists on their rotation 
amplitude or orientation. Initial models harboring rotating 
bulges would not only modify the inner parts of the rotation 
curves, but also the final mass distribution, due to a differ- 
ent spin-orbital coupling. New models would be required 
to quantify the strengths of these effects. 

4. The orbit. 

We use a single energy and inclination for the merger or- 
bits. While most aspects of the merger transformations 
should remain unchanged under variations in the orbital 
configurations, differences in the details are to be ex- 
pected. Sampling a complete set of preferential orbits 
according to the distributions extracted from cosmologi- 
cal simulati ons would allow us to refine our q uantitative 
_predictions fearitskv & Rix I ll997[ l 7aritskv et al...1997a.b^ 
Ikhochfar & Burkert 200^. 

5. Gas and star formation. 

Our models do not include gas hydrodynamics and star for- 
mation. One may argue that the gas must not affect the stel- 
lar dynamics, given that gas accounts for less than 10% of 
the baryonic mass of a typical spiral galaxy. Hence, the 
gas "is there for the ride". However, changes in the dy- 
namics could occur as gas concentrating in the galaxy nu- 
cleus renders the inner potential more centrally-symmetric, 
and modifies the allowed orbits for stars and gas. With or 
without such effects, the gas would be highly responsive 
to gravitational perturbations, and any ensuing star forma- 
tion would alter the stellar make-up of the final bulges and 
dis cs. Changes of M /L, perhaps by as much as a factor of 
~2 jPortinari et alJ2004) . would modify our predictions for 
global scahng relations. Tests performed by modifying the 
assumed M/L of bulges and discs, find that the decreasing 
trend of log(re,Bui)-AlBui (Fig.ll Ih) reverses sign when in- 
creasing the bulge's M/L by a factor of 2. 

6. The presence of bars and other galactic components. 

Our primary galaxy differs from the vast majority of disc 
galaxies by not hay ing a central bar {Eskridge et al. 2000; 
iKnapen et alJ2 00(7). Initial central components such as bars 
or lenses, would alter the galactic potencial in the center, 
and again, the mass distribution might be altered. In partic- 



ular, a central bar might increase the mixing of bulge and 
disc material above that given by the merger. Simulations 
with bars in the initial models would be required to quantify 
the effect. 

7. The reference epoch for our primary galaxy model 

Because we model our primary galaxy after the MW, one 
may argue that our experiments only test recent bulge 
growth; and, similarly, that the importance of satellite ac- 
cretion for the growth of present-day bulges would bet- 
ter be studied using galaxies modelled after those that ex- 
isted when the minor merger rate was relevant. A specific 
question is whether enough galaxies had well-developed 
discs for the processes modeled in this paper to be astro- 
physically relevant bulge-formation mechanisms. One dif- 
ficulty is to determine what this reference epoc h might be- 
According to CDM cosmolo gical simulations jBertschikl 
2004; Bertschik & Burkert 2004) as well as to observations 
lILe Fevre et aljil200Qi) . the minor- and major-merger prob- 
abilities both peak around z > 1.0; although the position 
of the peak is very dependent on t he environment and th e 
mass of object one is looking at ('Gottlober et al.ll2 00ll) . 
At this epoch, the Hubble sequence was in plac e, galaxy 
sizes were broadly similar to todav's Isimard et a]Lpi999l: 
iTruiillo & A^uerrill2004 iRavindranafli et alJl2004 . and, 
therefore, our modeling would be applicable. On the other 
hand, several lines of evidence point to 1 < z < 2 as a key 
epoch for the appearance of present-day early-type ga laxian 
components, i.e., e llipticals and bulges (Lilly et ^[l^sS 
lEhche-Moral et a l. 2006). As we move to higher redshifts, 
dominant morphologies depart more and more from mature 
Hubble types. Galaxies were smaller, gas-rich, less geomet- 
rically regular, and were fo rming st ars at higher rates than 
local galaxies dSteidel et al. 1996; Lowenth al et'al]ll997l: 
IPriver et al. 1998; Bunker et al. 2000; Truiifl o et alJ2004h 
Few of the bulges formed at these early cosmic epochs (e.g. 
iKoo et alJEoOS) may have been put in place through the 
processes modeled here. 

5. Discussion 

We have shown that the accretion of a satellite simultane- 
ously rises B/D and the bulge n index. More generally, ac- 
cretion causes a systematic evolution that makes the final 
structural parameters of the galaxy to be strongly correlated 
with the bulge luminosity. Given that the infall of small, col- 
lapsed baryonic systems is an inh erent ingredient of galaxy 
formation mode ls based on CDM ( Ka uffmann & Whitell 19931: 
iNavarro & Steinmetz..200Q,) . it is natural to ask what the ef- 
fect of many such accretion events might be. Considering that 
the effects of merging are probably cumulative (ABPOl), the 
lessons derived from the models presented here suggest that 
bulge growth and n increase would be progressive, leading to 
a continuous evolution of disc galaxies toward earlier Hubble 
types. 

KK04 list satellite accretion among the "environmen- 
tal, secular evolution" processes that drive galaxy evolution. 
Indeed, as long as satellites galaxies are available, satellite ac- 
cretion should be recurrent over long timescales. Even a single 
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accretion event may be termed a secular process given its dura- 
tion of well over the galaxy's orbital time (our merger models 
take ~2 Gyr to complete), and the low amplitude of the grav- 
itational perturbation by the satellite. But KK04 address the 
role of accretion in (pseudo)bulge growth by focussing on the 
ability of accretion to furbish gas for star formation. Our work 
suggests that accretion events may also contribute to bulge 
growth through an altogether different process, namely the in- 
ward transport of disc stars to the region of the bulge. 

The products of satellite accretion may even share some of 
the properties that are commonly used to defin e pseudobulges 
(as listed by, e.g., iKormendv & Fisher1l2005l) . and which are 
usually blamed on disc instabilities. Rapidly-rotating bulges 
are expected from prograde accretions, which deposit angular 
momentum in the bulge region and increase any pre-existing 
rotation (compare before- and after- primary bulge rotation 
curves in Figure|9}. Furthermore, inward piling-up of disc ma- 
terial in the region of the bulge should add rapidly-rotating 
material to the region of the photometric bulge. Inner spi- 
ral patterns are also expected from accretions, especially so 
in low-B/D systems in which gas from the primary or from 
the satellite settle into an inner disc. The limited resolution of 
our models prevents us from making quantitative statements. 
Measurement of a bulge rotation velocity not afifected by the 
inner disc is hard given the thickness of the disc. Measuring an 
ellipticity for the bulge is hard for the same reason. 

X-/peanut-shaped structures, another definitory feature for 
pseudo-bulg e s, are not seen in any of our models. However, 
Mihos etall (1 19951) show an accretion model which results 
in a strong X-shaped structure. Their primary galaxy, which 
lacks a bulge component, develops a strong bar which buckles 
into a peanut-shape. In our models, the presence of a massive 
bulge in the primary galaxy stabi lizes the disc against a strong 
bar response ( Sellwood & Moore 1 19991 bebattista et alJ2005[ 
GGB05). Peanut-shaped bulges may thus appear after the ac- 
cretion into low-B/D galaxies. 

The above arguments suggest that satellite accretion may 
play a role in setting those signatures commonly used to define 
pseudobulges. The role of the accretion is that of a trigger of 
the disc perturbation which leads to the infall of disc material 
to feed the bulge. 

In fact, satellite accretion events can account for several ob- 
servational facts that disc internal instabilities can not. Firstly, 
they can connect external and internal processes, as several 
works h ave addressed that bars are ea sily excited by tidal inter- 
actions llNoguchil[l987l: iGerin et alJ h99fl: Miwa & Noauchi 
And secondly, minor mergers are in agreement with the 
statistical association of young populatio ns wi t h mor pholog- 
ical disturbances reported by Kannappan et alJ ll2004 : while 
being, at the same time, also compatible with the low fraction 
of very peculiar/distorted shapes found at z <0.4 (<10%, see 
Ivan den Bergh et al.l2000l) . Satellite accretion events can be the 
key processes to reconcile the high merging rates of hierarchi- 
cal formation scenarios and the fact that many disc galaxies do 
not exhibit signs of major merger violence. 



6. Conclusions 

Our simulations of TF-scaled satellite accretion experiments 
show that: 

1. The TF-scaled satellite accretion drives pre-existing discs 
to a secular evolution towards earlier types (higher B/D and 
n). More massive satellites produce higher changes in n, 
B/D ratio, and disc thickening. 

2. For the BjD-Q.S primaries modeled here, BjD and n grow 
such that, roughly, n increases by a factor of ~ 2.3 when 
BjD doubles. Steeper growth of n is expected in primary 
galaxies with smaller initial bulges. 

3. The evolution of the bulge component is driven foremost 
by the injection of disc material to the center through tran- 
sitory non-axisymmetric distortions, and to a lesser extent 
by expansion of the pre-existing bulge and by deposition of 
satellite material. Vertical heating thickens the disc. 

4. The tidal radius of the satellite depends critically on the 
choice of TF exponent assumed for the scaling of galaxy 
densities. For the B/D-0.5 primaries modeled here, ajp - 
3 leads to complete disruption and to a torus of accreted 
material, while the satellites reach the center for axF = 4. 

5. The satellite disc deposits its mass onto the remaining pri- 
mary disc. Particles from the satellite bulge settle into a 
rather cold inner disc or ring, due to orbit circularization. 

6. The radial and vertical redistribution of material by heating 
produces a population mixture. Part of the material initially 
belonging to the primary disc is interpreted as bulge in the 
decompositions. The photometrically identified final bulge 
is partially supported by rotation, as it is observed in several 
galaxies. 

7. Final rotation curves exhibit symmetric shapes, and a gen- 
tle outward decline typical of early-type disc galaxies. The 
satellite material does not dominate the central parts of the 
final remnants in any case, and hence, central counterrota- 
tion is not visible in the total rotation curves of retrograde 
mergers. 

8. Accreted TF-scaled satellites cause systematic structural 
transformations in the primary galaxy, which evolves to- 
wards higher B/D, higher n, higher ctq, higher /iDisc and 
lower yUo,Disc, all following trends with the bulge magni- 
tude similar to the correlations observed by BGP04b. On 
the other hand, some observational trends that are observed 
to not correlate with the bulge magnitude do it strongly in 
our remnants, as the ratio between the central intensities 
of the bulge and the disc, and the bulge eff'ective surface 
density. We fail to reproduce the trends involving the bulge 
effective radius re.Bui- 

9. Secular evolution of galaxy bulges, currently being dis- 
cussed as the outcome of bar instabilities in the disc due 
to gas cooling, can be due to satellite accretion as well, 
through the inward transport of disc material to the region 
of the bulge. 
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